We present results from a nonlinear, 3D, time dependent numerical spectral model (NSM), which extends from the ground up into the thermosphere and incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). Our focal point is the mesosphere that is dominated by wave interactions. We discuss diurnal and semi-diurnal tides ill the present paper (Part 1) and planetary waves in the companion paper (Part 2). To provide an understanding of the seasonal variations of tides, in particular with regard to gravity wave processes, numerical experiments are performed that lead to the following conclusions: 1. The large semiannual variations in tile diurnal tide (DT), with peak amplitudes observed around equinox, are produced primarily by GW interactions that involve, in part, planetary waves. 2. The DT, like planetary waves, tends to be amplified by GW momentum deposition, which reduces also the vertical wavelength. 3.Variations in eddy viscosity associated with GW interactions tend to peak in late spring and early fall and call also influence the DT. 4. The semidiurnal semidiurnal tide (SDT), and its phase in particular, is strongly influenced by the mean zonal circulation. 5. The SDT, individually, is amplified by GW's. But the DT filters out GW's such that the wave interaction effectively reduces the amplitude of the SDT, effectively producing a strong nonlinear interaction between the DT and SDT. 6.) Planetary waves generated internally by baroclinic instability and GW interaction produce large amplitude modulations of the DT and SDT.